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Abstract 

The '^^S(p, 7)^^01 reaction is expected to provide the dominant break-out path from the SiP 
cycle in novae and is important for understanding enrichments of sulfur observed in some nova 
ejecta. We studied the ^^S('^He,t)'^^Cl charge-exchange reaction to determine properties of proton- 
unbound levels in '^^Cl that have previously contributed significant uncertainties to the '^^S(p,7)^^Cl 
reaction rate. Measured triton magnetic rigidities were used to determine excitation energies in 
^^Cl. Proton-branching ratios were obtained by detecting decay protons from unbound '^^Cl states 
in coincidence with tritons. An improved ^^S{p,^)^'^Cl reaction rate was calculated including robust 
statistical and systematic uncertainties. 



* E-mail: lmatos@lsu.edul 
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I. INTRODUCTION 



A classical nova is a cataclysmic nuclear explosion on the surface of a white dwarf star 
resulting from the accretion of hydrogen-rich gas from a companion star. A sequence of 
nuclear reactions produces a sudden luminosity increase by up to a factor of a million and 
ejects matter from the white dwarf. The time scales of explosive hydrogen burning processes 
are influenced by the duration of reaction cycles closed by {p,a) reactions, with break out 
via (p,7) reactions competing with /3-decays. The SiP cycle is one such cycle, which is of 
particular interest for understanding novae like Nova Her 1991 that are observed to exhibit 
high sulfur abundances compared to solar values 

HQ. 



In the SiP cycle the S(p,7) CI reaction is believed to be the dominant break-out re- 
action 3|. The rate of the '^^8(^,7)'^^ CI reaction at nova temperatures is dominated by 
resonances corresponding to states in the compound nucleus ^^Cl. The ^^S(p,7)^^Cl reaction 
rate was previously calculated based on some measured resonance properties and estimates 
of others based upon the mirror nucleus ^^P {4]. Uncertainties in the rate were provided in 
a subsequent reanalysis |5|. The ^^Cl excitation energies near the proton threshold of Qp 
1581.3(6) keV 6] used by Ref. 
the 2 Sigma level (by 10-20 keV) 



5 1 were based on earlier measurements that disagreed at 
Q], and the resonance strengths were only constrained 



based upon properties of states in the mirror. A recent study of 5 states below 2.3 MeV 



n 



via the 32s(3He,t)32ci reaction [8] is in agreement with the excitation energies reported by 
Ref. j^. However, significant uncertainties remain regarding resonance strengths and size 
of systematic uncertainties in the resonance energies. 



We have studied proton-unbound states in ^^Cl using the ^^S('^He,t)'^^Cl charge-exchange 
reaction. Excitation energies and proton-branching ratios for states of astrophysical interest 
for the ^^S(p,7)^^Cl reaction rate were determined. In the following sections we describe 
the experiment and new results for states in ^^Cl, including the observation of a predicted, 
but previously unobserved, level. We then present calculations of a new recommended 
^^S(p,7)'^^Cl reaction rate based on these new results including statistical and systematic 
uncertainties. 
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FIG. 1. Spectra showing position of tritons at the focal plane measured by the ionization chamber 

with the Enge spectrograph set at 3°. Three spectra show results from measurements with three 
targets, 240 //g/cm^ ZnS, 350 jig/cm^ ZnS and a 300 /xg/cm^ Si target. Peaks are identified with 
the final level in the nucleus produced. 



II. EXPERIMENT 

We used the charge-exchange reaction ^^S(^He,i)^^Cl to populate states in ^^Cl. A 30- 
MeV ^He^+ beam from the Extended Stretched TransUranium (ESTU) tandem Van de 
Graaff accelerator at the Wright Nuclear Structure Laboratory (WNSL) at Yale University 
bombarded ZnS targets. Targets with thicknesses of 240 /xg/cm^ and 350 /ig/cm^ both on 5 
//g/cm^ carbon substrates were produced via evaporation at Oak Ridge National Laboratory 
and used in the experiment at the WNSL. Their thicknesses were determined with about 
10% uncertainty via energy loss measurements of a particles from a ^^^Am source. A 150 
jig/cm^ CdS target on a 20 jig/cw? substrate was used at one angle (5°) but had worse energy 
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resolution and lower counting rates and was therefore not used at other angles. Additional 
data were taken with a 300 yUg/cm^ natural Si target for calibration and a 900 /xg/cm^ Zn 
target for background analysis. 

Reaction products were separated and analyzed with the Enge split-pole spectrograph at 
the WNSL. The spectrograph separates particles according to their magnetic rigidities, Bp, 
so that discrete positions at the focal plane correspond to discrete momenta. The positions of 
detected particles at an Enge focal plane were determined by a position-sensitive ionization 
drift chamber filled with 150 Torr of isobutane gas 91]. As the ions pass through the gas, the 
amount of charge collected by the cathode determines the energy lost in the detector. The 
residual energy was measured by a thick plastic scintillator located behind the ionization 
chamber. Measurements were conducted at spectrograph (laboratory) angles of 3°, 5°, 10° 
and 20°. At the 3° setting the protons emitted from excited states in '^^Cl were detected in 
coincidence with outgoing tritons by the Yale Lamp Shade Array (YLSA) |l^ consisting of 
four 16-strip silicon detectors arranged in a lamp-shade configuration. The YLSA detectors 
covered an angular range of 6iab = 131° to Oiab = 166° and were calibrated with a particles 
from a ^^^Am source. 



III. LEVEL ENERGIES 

Tritons were identified at the focal plane of the spectrograph using the Eres vs. AE 
relationship from the ionization chamber and the scintillator. Focal-plane position spectra 
gated on tritons are shown in Fig. [H and states in '^^Cl populated via the ^^S('^He,)f:)'^^Cl 
charge-exchange reaction are labeled. Tritons corresponding to states in ^^F resulting from 
the i60(3He,t)^*'F reaction due to oxygen contamination in the targets are also seen in Fig. [H 

Positions at the focal plane were calibrated using known states in populated from the 
^^Si(3He,t)2^P reaction on a silicon target. The 424 keV ^^F state populated from oxygen 
contaminants in the targets, and the ground and 90 keV first excited states in ^^Cl were 
also used to calibrate the position spectra, providing 14 calibration points spanning the 
region of interest. Centroids were obtained by fitting a Gaussian function to each peak with 
a linear background, and the centroids of the calibration peaks were fit by their magnetic 
rigidities. Both 2"'' and 3'"'^ order polynomial functions were fit. Excitation energies for 
all states of astrophysical interest agreed within 1 keV from both fits, but the 3'''^ order fit 
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function was adopted as it provided a better fit to the lowest and fiighest energy calibration 



peaks. The magnetic rigidities of the recoiling tritons were ca 
kinematics with the newest atomic mass values implemented 11 



culated from the reaction 



12[, including new values 



for and ^201 obtained recently, with mass excess of -7147.5(12) keV, -13333.8(12) 

were made for the energy loss of the 



incident '^He ions and recoiling tritons in the target using the energy loss code STOPIT 



assuming the reaction happened in the center of the target. The fit function was then used 
to obtain magnetic rigidities from the centroids for the peaks of interest. The energies of 
the levels were calculated from reaction kinematics taking energy losses into account. The 
statistical uncertainty in the energy of each state, AE'j, was calculated as a combination 
of the uncertainties originating from the centroid determination and the fitting function 
uncertainties estimated from covariance analysis. The energy resolution for each triton peak 
was typically about 40 keV (FWHM), and no peaks wider than the instrumental resolution 
were observed. 

Excitation energies were obtained for each level at several angles and with the different 
targets. The final weighted average for each state was calculated from these energies. The 
uncertainty in the excitation energy was assigned as the larger of the average uncertainty 



(1) 



and the scatter uncertainty 



AE 
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Systematic uncertainties were estimated using a Monte Carlo simulation, where the target 
thicknesses and mass values were randomly varied with Gaussian distribution probabilities. 
Mass and thickness uncertainties were used as standard deviations, a. The systematic error 
was determined as the standard deviation of the Gaussian fit for the level energy distributions 
obtained in the Monte Carlo simulation. The systematic uncertainty was found to be 4 keV 
for all states, except for the 462 keV and 1167 keV states, where the systematic uncertainty 
was determined to be 2 keV and 3 keV, respectively. 
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FIG. 2. (Color online) Focal-plane position spectrum with the Enge spectrograph set at 10° showing 
the newly discovered 2611 keV state in ^^Cl produced in the 350 /xg/cm^ ZnS target. 



The excitation energies determined from this measurement with statistical uncertainties 
are given in Tabled They are com^red with other measurements 0, 7, 8] and the evaluation 



by 



measurements were 



ISH . The evaluation favored 3| over Ml, where results from both 
available and generally adds 11 keV to |7| if only this one is available. The level energies 
determined in ^^Ar(/3+p) experiments include 461.1(1) keV, 1168.5(2) keV, and 4076(10) keV 
from 13] and 4072(9) keV from [l7], all agreeing well with our results. A recent measurement 
also using the ^^S(^He,t)^^Cl reaction provides energies for five levels between 1.3 and 2.3 
MeV that are systematically higher than the current work by about 4 keV 



We have also discovered an excited state in ^^Cl at 2611(5) keV, shown in Fig. O The 
measured excitation energy for the new state is in good agreement with a l"*" state predicted 
at 2574(50) keV based on the mirror nucleus ^^P and estimated using the IMME equation 4|. 



IV. PROTON UNBOUND STATES 



The YLSA silicon-detector array was installed covering backward angles at the target 
position to measure decays of proton-unbound states in ^^Cl in coincidence with tritons 
detected with the spectrograph at 9iab = 3°. Figure [3] shows the energy of the protons 
detected by YLSA versus the position (corresponding to momentum) of the tritons at the 
focal plane. A gate is shown around events corresponding to the proton decay of ^"^Cl* to 
the ground state of ^^S (J'" = 1/2+). Proton decay to excited states of ^^S was not possible 
for ^^Cl excited states below 2.8 MeV due to the '^^Cl proton-separation energy of Qp = 
1581.3(6) keV and the first excited state in ^^S being at 1248.9(2) keV Q. The events 
below the gated region in Fig. [3] corresponds to the proton decay of ^^F*, while events above 
are caused by leakage of deuterons into the triton window in the particle identification cut. 
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FIG. 3. (Color online) Triton energy - proton energy coincidence spectrum. The gated band 
corresponds to the ^^S(^He,i)^^Cl*(p)^^S reaction. 
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FIG. 4. (Color online) Triton-proton angular correlation probabilities from the '^^S(^He,t)^^Cl(p)'^^S 
reaction at 30 MeV for the various states in '^^Cl as listed in Table HI Red squares are the exper- 
imentally determined values with uncertainties, the black lines are the fits with even Legendre- 
polynomial terms. The polynomial orders are listed in Table HIl Although the 2"'^ order Legendre 
polynomial (dashed line) for the 2859 keV and 3659 keV states describes the data well, we assume 
the 4*'' order Legendre polynomial to be correct (solid line). 



The angular probability distributions of the emitted protons for the ^^Cl states between 
2.1 and 3.9 MeV are shown in Fig. |H Two neighboring YLSA strips are coupled together to 
reduce scatter due to the low statistics. The values come from the ratios of the number of 
events in the t-p coincidence peaks to the total number of tritons populating a given state 
in ^^Cl, both with background subtracted. These numbers were then divided by the YLSA 



io| 



efficiency estimated via the use of a Monte Carlo model 

Angular correlations of isolated nuclear levels can be described by a linear combination 
of even terms of the Legendre polynomials, P2k [cos{9c.m.)], of center-mass angle Oc.m. up to 
two times the proton orbital angular momentum i, i.e. 

I 

W{9) = Y.AkP2k[cos{ec.m.)], (3) 

k=0 

which is symmetric around 6 cm. = 90°. The minimum order of the Legendre polynomial 



needed to fit the data is defined by statistical significance testing with a p- value 



19| required 



to be > 0.05, corresponding to < 14.07, 12.59, and 11.07 for 7, 6, and 5 free parameters, 
respectively. This function was constrained to be positive at each point, and the integral 
over the full solid angle must be < 1. The function was then integrated over the full solid 
angle to obtain a total proton-branching ratio, bp, as 

bp= r 2-Ksin{e)W{e)d9. (4) 
Je=o 

The results are shown in Table ITTl with Pfu being the order of the Legendre polynomial used 
in the fit. The quoted proton branching-ratio uncertainties come from the uncertainties 
in the fit parameters. 

The proton-branching ratio for the first excited state above the proton separation energy. 
Ex = 1734.2 keV, is expected to be small. However, even if it was a significant branch, 
the proton energy resulting from the Ej. = 153 keV resonance is below the YLSA detector 
threshold, and would not be observed. The angular distributions of states with excitation 
energies of 2128 keV, 2203 keV, 2279 keV, 2611 keV, 2677 keV, and 2932 keV are well fit 
with just a constant function, requiring only the first term in Eq. [3l 

The new state at 2611 keV presents an interesting case. Fits are shown in Fig. H] using 
only an isotropic term (solid line) and the 4*'^ order Legendre polynomial (dashed line). 
The angular distribution for this state is better fit with the 4*^* order Legendre polynomial, 



but the isotropic fit fulfills the p- value test [l9| at the 95% confidence level and cannot be 



ruled out. As the excitation energy agrees well with that expected for a state corresponding 
to a 1+ state in the mirror [4], we adopt the 1+ assignment and use the isotropic angular 
distribution required from a £=0 proton orbital angular momentum. 

For the 2859 keV state, the 2"^^ order Legendre polynomial fulfills the p-value test giving 
a total proton-branching ratio of 75±5 %. However, a fit with the 4*^* order Legendre 
polynomial (dashed line in Fig. H]) differs from the previous one only outside the area covered 
with our data points, giving a total proton-branching ratio >95%. The neutron spectroscopic 
factor for the mirror state in ^^P has been measured to be 0.03 [22:] and 0.008 [23,]. While 
there is a discrepancy between the measurements, even the lower value implies an expected 
proton width for the 2859 keV level that would be about 3 orders of magnitude larger than 
the expected gamma width. Therefore, we adopt the result from the 4^^ order Legendre 
polynomial fit. 

The situation is similar for the 3695 keV state. We find a total proton-branching ratio 

10 



of 59±4 % coming from the fit with the 2"'^ order Legendre polynomial and a branching 
ratio of >96% using a 4*'' order Legendre polynomial. While the mirror assignment is not as 
clear for the 3695 keV state, the most likely candidate, the 3880 keV (2+), has a measured 



spectroscopic factor of 0.028 22|, in agreement with ^0.03 predicted for the '^^Cl state by 
shell model calculations {24]. This is a factor of 15 more than the spectroscopic factor 
required for the branching ratio of 59% and we therefore adopt >96% from the 4*^^ order 
Legendre polynomial fit. 

For the 3055, 3164, 3281 and 3875 keV states, a 4*'* order Legendre polynomial fit is 
required, and the resulting branching ratio is consistent with bp = 100%. The lower limit 
for these states is then statistically estimated based on the number of events. 

The obtained values for minimum proton orbital angular momenta are in a good agree- 
ment with assumed spins and parities. The spin-parity assignments for states with < 
3.5 MeV are taken from j4|. The spins of the 3695 and 3875 keV states were tentatively 
assigned based on the mirror symmetry corresponding to the 3880.3 keV 2"*" and 3989.8 keV 



(3^) states in ^^P ISj, as the 3796.1 keV (1^) in ^^P is assumed to be the mirror state to 



the 3767 keV state in ^^Cl that is known to be = 1+ [16|. 



V. 3is(p,^)32ci REACTION RATE 

At nova temperatures, T ;^0. 1-0.3 GK, the states just above the proton-separation energy 
[Qp = 1581.3 keV) dominate the ^^S(p,7)'^^Cl reaction rate. As the resonances are generally 
narrow and well separated, the resonant component of the reaction rate (in cm^mol~^s~^) 
can be approximated by 

Na{(tv) = 1.54 X 10"(/iT9)"^/2 

X 5](a;7),exp(-11.605£;,/T9), (5) 

r 

where Tg is the temperature in GK, is the energy of the '^^Cl resonance in MeV, /i is the 
reduced mass in atomic mass units, and {uj'j)r is the resonance strength in MeV, given by 
the spin of the resonance, Jr and its partial (rp,r^) and total (F) widths as 

= l^I^IZl. (6) 

The resonance reaction rate depends exponentially on the resonance energies, Er, and lin- 
early on the partial widths through the resonance strengths, though the proton partial width, 
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Tp, also has an exponential dependence on energy through the penetrability. Therefore, our 
improved measurement of the resonance energies and the proton-branching ratios, corre- 
sponding to Fp/r = rp/(rp + r^), has a direct impact on the uncertainty in the ^^S(p,7)^^Cl 
reaction rate. 

We have calculated the gamma widths for the states in ^^Cl using the known half-lives 
of mirror ^^P states, T1/2, as well as 7-branching ratios, b^, and energies, E^., of the corre- 



sponding transitions from the states [20|] . Assuming that the reduced transition probabilities, 
B[Ei) and B{Mi), are the same for both mirror nuclei, one can calculate the 7 width of a 
state in the mirror nucleus as a sum through all possible final state transitions: 

P /32pA ^ V ^t'-^^H'^Cl) &,;.ln(2) 

where A is electric or magnetic multipolarity. This follows a similar prescription as was 
used in 0, 0]. The lowest possible multipolarities were assumed. In the case of M1/E2 



transitions, studies of the mirror nucleus 20|] showed that Ml transitions mostly dominate, 
and thus Ml transitions were adopted in the present reaction rate calculations. 

For the excited states between 2.1 MeV and 2.3 MeV, where the proton-branching ratio 
was determined to be finite but less than 100%, we calculated the proton widths, Tp, directly 
from the gamma widths with our measured proton-branching ratios. For the higher energies, 
the resonance strength becomes insensitive to the proton width as Tp ^ and FpF^/F ~ 
F,. 

The 1734 keV state, corresponding to the 153 keV resonance, is the one state where 
the proton decay width is important, but where no information was extracted from our 
measurement. In this case we calculated the proton width using the prescription also followed 



infl. 



Tp = 2^P,C'Sp9%, (8) 



where /i is the reduced mass and Oc = 5.6 fm is the channel radius. The single-particle 
reduced width, = 0.32, was derived from the parameterization for nuclei with the mass 
number A = 12-50 and bombar ding energies E < 1000 keV based on optical-model computa- 
tions and i?-matrix expressions |2l|. The penetrability was calculated to be Pc = 2.9 x 10~^^. 
The spectroscopic factor, C^Sp, was obtained from the reaction studies with the mirror nu- 
cleus '^^P produced via a neutron transfer, ^^P((i,p)^^P, at deuteron energies of 10 MeV 
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22| and 20 MeV 23|]. The spectroscopic factors reported in these measurements are dis- 

5|. However, 
23| using 



crepant, 0.011 j22| and 0.0054 j23|, and the average value was adopted in 
we conducted a reanalysis of the experimental cross section data from 



and 



the FRESCO code 25|. We find differential cross-sections from both experiments to be best 
fit with a spectroscopic factor of 0.011, in agreement with 22|. Thus, we adopt the higher 
value for the spectroscopic factor and from this calculate the proton width of the 153(5) keV 
resonance to be 2.8x10^*^ meV. 

The calculated proton and gamma widths, as well as resonance strengths calculated using 
Eq. [6l are listed in Table IIIII The resonance energies determined in this work have been 
used for the calculations, except for the 3767 keV state measured in the ^^Ar /3-decay studies 

n n 

|16l | and the 3397 keV state predicted by based on the mirror symmetry. 

To cover temperatures below the regions dominated by the resonances, we have adopted 
the direct capture parameterization from j4|. The recommended total ^^S(p,7)^^Cl reaction 
rate for the stellar temperature range T = 0.01-10 GK is given in Table HVl The individual 
contributions of direct and resonant capture are illustrated in Fig. [51 The rate has been 



parameterized in the Reaclib format 
of 7 parameters in the form 

NA{crv) = exp 



26| as the sum of three exponentials, each with a set 



+04X9 + a^Tg^^ + aQ X ln(Tc 



(9) 

where Tg is the temperature in GK, using the tools in [27[. The fit agrees with the data 
within 1.5%. The new Reaclib parameters are listed in Table IVl 

Uncertainties in the reaction rate have been estimated using a Monte Carlo technique, 
as a combination of normal and log-normal distributions of uncertainties complicate the 
analysis. Uncertainties in the resonance energies contribute to the log-normal distribution, 
as the energy is in the exponential of Eq. [51 The gamma width contributes to the normal 
distribution, as its uncertainty originated from the half-life uncertainty, mirror-symmetry 
assumption (that we have estimated to be 20% based on the mirror states in the neighboring 
nuclei) and the uncertainty of the proton-branching ratio. The uncertainty in the proton 
width will have a pure normal distribution, if extracted from the proton-branching ratio. 
For the 153 keV resonance, the exponential dependence of the penetrability on the energy 
contributes to the log-normal distribution, but a normally distributed contribution also 
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FIG. 5. (Color online) ^^S(p,7)^^Cl reaction rate as a function of the stellar temperature T. 
Resonances with a contribution of at least 10% contribution are shown. 

originates in the spectroscopic factor uncertainty that we have estimated to be 30%, the 
spectroscopic factor being different for both mirror nuclei due to the effect of Coulomb and 
other isospin-nonconserving interactions (~10%), a small uncertainty of the wavelength in 
the penetrability (~2%), and the reduced width uncertainty that we have estimated to be 
20%. For the 3695 keV and 3875 keV states, the uncertainty in the spins is considered. 
Estimated uncertainties for the resonance strengths are shown in Table IIIII Direct capture 
uncertainties were taken directly from j^l. 

In the simulation, the values of the energies and resonance strengths were varied randomly 
as a Gaussian distribution. Correlations, when dependent on the same parameters, were 
taken into account. The resulting reaction rate distributions (see Fig. |6]) have various shapes 
for different temperatures, including a nearly pure log-normal distribution at T = 0.05 GK 
and normal distribution at T = 5 GK. To give final uncertainties that would correspond 
to the standard deviation a in both distributions, we found a lower limit as a value with a 
percentile of 15.9 and an upper limit with a percentile of 84.1, covering the 68.2% confidence 
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level. The results from the Monte Carlo simulation are listed as the low and high rates in 
Table Ml 



VI. DISCUSSION 

The total ^^S(p,7)^^Cl reaction rate and individual contributions based upon this work 
are illustrated in Fig. |5l Direct capture dominates the reaction rate up to T ~ 0.03 GK. 
The 153 keV, 546 keV, 622 keV, and 1096 keV resonances, corresponding to 1734 keV, 2128 
keV, 2203 keV, and 2677 keV levels in ^"^Cl, dominate the rate over nearly all temperatures. 
The 697 keV resonance contributes more than 10% at nova temperatures, and the 1030 keV 
resonance must be considered at X-ray burst temperatures, T ^ 2 GK. The 2186 keV, and 
2294 keV resonances do not contribute except at very high temperatures, T > 5 GK, and 
the 2859 keV and 3695 keV levels (which had some ambiguity as to the shape of the proton 
angular distribution) make a negligible contribution to the reaction rate. 

The ^^S(p,7)^^Cl reaction rate with uncertainties was recently calculated by using the 
previous work of [4] with the evaluated level energies from |l^. The ratio of the rate from 
Ref. [5i] to our rate is shown in Fig. [71 The uncertainties in both rates are illustrated by 
the hashed regions. The agreement for T < 0.03 GK is expected, as the direct capture rate 
was calculated based on the same parameters and the lower and upper limits we take 
from [5]. 

Over much of the range of nova temperatures our recommended rate is significantly 
greater than even the "high rate" recommended in Ref. |5|. This arises from the contribution 
of individual resonances. In Fig. [8] the individual resonance reaction rates from Ref. [Sj] are 
compared to our results. Our higher reaction rate at most nova temperatures arises from 
the fact that the resonance energies adopted in Ref. J] (and derived from 15|) are greater 
than our energies by 6-15 keV (9 keV on average). The excitation energies adopted by 
[isl ] and 4] primarily reflect a weighted average of 3] and j3] after the results of [3] were 
shifted to match the 1168 keV excitation energy. The 3-5 keV uncertainties in the adopted 
excitation energies do not properly reflect systematic uncertainties, such as uncertainties in 
the calibration used or the discrepancies between measurements. 

Our rate also differs from that of [s^ due to improved values for resonance strengths. It 
should be noted that while we adopt a resonance strength for the 153 keV resonance that 
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is nearly the same as {4, 5], this arises from two significant changes that largely cancel each 
other. Our lower resonance energy results in a significantly smaller penetrability, but we 
recommend a significantly greater spectroscopic factor based on our reanalysis of the transfer 



data of 



22 



23|. Our reaction rate near peak nova temperatures, T ^ 0.3 GK, becomes 



smaller than that of js] due to our improved values for the proton branching ratios of the 
546 and 622 keV resonances. We find r.„ for the 546 keV resonance to be about 30% smaller 



than estimated by j4| and Tp for the 622 keV resonance to be about 2 orders of magnitude 
smaller. The smaller proton widths result in smaller resonance strengths and reaction rates, 
though the magnitude of the effect is mitigated by the lower value of the resonance energies 
for these states that makes for a smaller decrease in the reaction rate than would otherwise 
be. 

At higher temperatures (above about 2 GK) our rate increases in comparison to jj] due 
to our inclusion of resonances above 2 MeV that were not previously considered. 

In summary, we have significantly improved the resonance energies and resonance 
strengths for some of the most important resonances in the ^^8(^,7)'^^ CI reaction. An 
important aspect of the current work is that we have given careful consideration to uncer- 
tainties, including systematic uncertainties in the level energies, states used for calibration, 
and target thickness effects. The largest uncertainties in the reaction rate at nova temper- 
atures arise from the systematic uncertainty in the resonance energies and the resonance 
strength of the 153 keV resonance. Our excitation energy for the 1734.2(14) keV state 
(corresponding to the 153 keV resonance) differs from the value of 1736.7(6) reported by 8| 
(which used a slightly different set of calibration reactions) by 2.5 keV or slightly more than 
la. However, we estimate the systematic uncertainty in the resonance energies to be 4 keV, 
which is in agreement with the fact that the excitation energies for all levels reported by js] 
are higher on average than this work by about 4 keV. Additional experimental information 
leading to an improvement in the resonance energies would therefore be valuable. As the 
states most important for novae have substantial branches for gamma decay, an accurate 
measurement of gamma-ray energies using a complementary approach, for example as in 
would be particularly helpful in reducing the systematic uncertainties that arise largely 



from Q-value uncertainties in reaction studies like this one. A direct measurement of the 
resonance strength of the 153 keV resonance (or corresponding proton width) is also desired, 
but would be experimentally challenging. 
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TABLE I. Excitation energies in CI measured in this work compared to other measm'ements and 
evaluations. Ah hsted uncertainties are statistical. Systematic uncertainties have been estimated 
in this work as 4 keV for all states, except for the 462 keV state with the systematic uncertainty 
of 2 keV and the 1167 keV state with the systematic uncertainty of 3 keV. 



this 


Jeanperrin Vouzoukas 


Endt 


Wrede 


work 


m 


[3] 


[15j 


[8] 


[keV] 


[keV] 


[keV] 


[keV] 


[keV] 


462.0(10) 


447(7) 




461.1(1) 




1167.0(21) 


1157(5) 




1168.5(2) 




1327.4(29) 


1326(5) 


1329(3) 


1331(3) 


1331.2(5) 


1734.2(14) 


1719(4) 


1735(3) 


1733(2) 


1736.7(6) 


2127.5(19) 


2122(5) 


2129(3) 


2130(3) 


2131.1(4) 


2203.1(28) 


2193(7) 


2213(3) 


2212(3) 


2209.5(5) 


2278.6(25) 


2270(5) 


2281(3) 


2281(3) 


2283.5(5) 


2610.9(30) 










2677.0(16) 


2665(10) 




2676(10) 




2859.2(14) 


2858(5) 




2869(5) 




2931.5(17) 


2941(5) 




2952(5) 




3054.7(14) 


3056(5) 




3067(5) 




3163.9(11) 


3166(5) 




3177(5) 




3280.8(23) 


3290(10) 




3301(10) 




3695.0(9) 


3692(7) 




3703(7) 




3874.8(17) 


3883(5) 




3894(5) 




3999.5(12) 


4002(6) 




4013(6) 




4073.6(11) 


4080(7) 




4074(7) 




4349.9(23) 


4356(7) 




4367(7) 




4577.1(30) 


4590(8) 




4601(8) 
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TABLE II. Proton-branching ratios for states in '^^Cl. Spin-parities assigned based on J], except 
for 3695 keV and 3875 keV states that we assigned based on the mirror symmetry. 



[keV] J" 


Pf^t 


hp [%] 


2128 


3+ 





7 ± 4 


2203 


1+ 





54 ± 7 


2279 


2+ 





66 ± 13 


2611 


1+ 





>62f. 


2677 


2+ 





>78^ 


2859 


3+ 


4 


>95 


2932 


2- 





>88 


3055 


4- 


4 


>97 


3164 


3" 


4 


>96 


3281 


2+ 


4 


>88 


3695 


2+ 


4 


>96 


3875 


3+ 


4 


>94 



b 



Fit result is (95±32)%. 
Fit result is (94±16)%. 
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TABLE III. Properties of proton-unbound states in ^^Cl and corresponding resonances relevant 
for the '^^S(p,7)^^Cl reaction. Energies (including systematic and statistical uncertainties) are the 
result from this work, except for 3767(10) keV from [16] and 3397(50) keV predicted by [J]. The 
spin-parity assignments from were used for the states up to 3.5 MeV, from jl^ for the 3767 keV 
state, for the 3695 and 3875 keV states were tentatively assigned based on the mirror symmetry. 



Ex 


Ef 




r. 






cj(u;7) 


[keV] 


[keV] 




[meV] 


[meV] 


[meV] 


[meV] 


1734 


153(5) 


3+ 


1.0 


2.8x10"^ 


4.9xl0~s 


1.0x10-®^ 


2128 


546(5) 


3+ 


7.9 


0.59 


0.96 


0.61 


2203 


622(5) 


1+ 


15.5 


18.2 


6.3 


3.6 


2279 


697(5) 


2+ 


3.1 


6.0 


2.54 


0.76 


2611 


1030(5) 


1+ 


20.2 




14.4 


7.5 


2677 


1096(5) 


2+ 


57.9 




68 


30 


2859 


1278(4) 


3+ 


5.4 




9.5 


3.7 


2932 


1350(5) 


2- 


2.3 




2.84 


0.76 


3055 


1473(5) 


4~ 


0.8 




1.81 


0.42 


3164 


1583(4) 


3- 


2 




3.51 


0.96 


3281 


1700(5) 


2+ 


15 




18.3 


8.5 


3397 


1816(50) 4+ 


1.8 




4.1 


1.4 


3695 


2114(4) 


2+ 


28 




35 


20 


3767 


2186(10) 1+ 


110 




83 


34 


3875 


2294(4) 


3+ 


59 




104 


65 



The uncertainty distribution of the 153 keV resonance strength includes an additional log-normal 
component with aiog{uj^) = 0.58. 
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TABLE IV. Recommended stellar reaction rates as a function of the temperature T for the reaction 
^^S(p,7)^^Cl. Lower and upper limits cover the 68.2% confidence level. 



Tempe- 


Recommended 


Low 


High 


rature 


rate 


rate 


rate 


T 


Na{(tv) 


Na{(tv) 


Na{(tv) 


[GK] 


[cm^mol~^s~"^] 


[cm^mol-is-i] 


[cm3mol-is-i] 


0.01 


3.68x10-^"^ 


2.33x10-^^ 


5.12x10-^^ 


0.015 


1.84x10-3^ 


1.17x10-37 


2.56x10-37 


0.02 


3.08x10-33 


1.95x10-33 


4.28x10-33 


0.03 


6.33x10-28 


4.36x10-28 


1.02x10-27 


0.04 


5.39x10-23 


1.39x10-23 


2.12x10-22 


0.05 


2.69x10-19 


8.50x10-20 


8.10x10-19 


0.06 


7.56x10-1^ 


2.86x10-17 


1.91x10-16 


0.07 


4.10x10-15 


1.76x10-15 


9.13x10-15 


0.08 


7.98x10-1^ 


3.77x10-1^ 


1.62x10-13 


0.09 


7.86x10-13 


4.00x10-13 


1.48x10-12 


0.1 


4.82x10-12 


2.59x10-12 


8.61x10-12 


0.15 


9.73x10-1° 


6.17x10-1° 


1.48x10-9 


0.2 


1.22x10-8 


8.42x10-9 


1.71x10-8 


0.3 


9.86x10-^ 


5.70x10-7 


1.45x10-6 


0.4 


1.41x10-^ 


7.90x10-5 


2.05x10-^ 


0.5 


2.99x10-3 


1.73x10-3 


4.26x10-3 


0.6 


2.27x10-2 


1.33x10-2 


3.20x10-2 


0.7 


9.47x10-2 


5.61x10-2 


1.33x10-1 


0.8 


2.73x10-1 


1.63x10-1 


3.81x10-1 


0.9 


6.17x10-1 


3.70x10-1 


8.59x10-1 


1 


1.18x10° 


7.16x10-1 


1.63x10° 


1.5 


8.39x10° 


5.59x10° 


l.llxlQi 


2 


2.37xl0i 


1.69xl0i 


S.OlxlQi 


3 


6.99x101 


5.25x101 


8.71x101 


4 


1.20x102 


9.08x101 


1.48x102 


5 


1.63x102 


1.25x102 


2.00x102 


6 


1.97xl02 


1.53xl02 


2.41 xl02 


7 


2.23xl02 


1.75xl02 


2.70xl02 


8 


2.41x102 


1.91xl02 


2.91x102 


9 


2.53x102 


2.02x102 


3.06x102 


10 


2.61x102 


222.08x102 


3.14x102 



TABLE V. Recommended Reaclib parameters for the '^^8(^,7)^^ CI reaction rate within T = 0.01-10 

GK. 

set ao ai a2 as a,^ as a^ 

1 297.13 -2.6702 114.07 -485.67 63.95 -5.699 140.32 

2 -35.362 4.1263 -324.5 373.08 -18.087 0.9159 -205.99 

3 1315.5 -1.8787 330.21 -1911.9 261.17 -25.381 511.02 
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FIG. 6. (Color online) Total reaction rate probability density functions as a result of the Monte 
Carlo simulation of the input-uncertainty propagation for various temperatures. The left figures 
show the distributions in the logarithmic scale normalized to the recommended value, the right 
figures show the same distributions in logarithmic scale without normalization. 
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FIG. 7. (Color online) The '^^8(^,7)^^ CI reaction rate with uncertainties calculated by [Sj (red), 
with only statistical uncertainties considered, is compared to the results from this work (blue), 
with both statistical and systematic uncertainties considered. Values are shown normalized to the 
recommended rate from this work. 
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FIG. 8. (Color online) The ratio of the individual resonance reaction rates from Ref. Q to rates 
from this work as a function of stellar temperature for the most important resonances contributing 
to the 3iS(p,7)32Cl rate. 
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